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Second electron affinities of Aujv and Agjv clusters and 
the dissociation energies for fission of the Au^~ and Ag^~ 
dianions are calculated using the finite-temperature shell- 
correction method and allowing for triaxial deformations. Di- 
anionic clusters with N > 2 are found to be energetically 
stable against fission, leaving electron autodetachment as the 
dominant decay process. The second electron affinities exhibit 
pronounced shell effects in excellent agreement with measured 
abundance spectra for Au^~ (N < 30), with appearance sizes 
n 2 a ~{Au)= 12 and n^(Ag)= 24. 
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Unlike the case of multiply charged cationic species, 
the production and observation of gas-phase doubly an- 
ionic aggregates had remained for many years a chal- 
lenging experimental goal. However, with the availabil- 
ity of large carbon clusters (which can easily accomo- 
date the repulsion between the two excess electrons) this 
state of affairs changed, including observation of doubly 
negative fullerenes,El C 60 ~ , and fullerene derivatives^ as 
well as a recent measurement oLthe photoelectron spec- 
trum of the citric acid dianionE Moreover, such obser- 
vations are not limited to carbon based aggregates and 
organic molecules, with a first observation of doubly 
ionic metal clusters (specifically gold clusters) reported 
most recently. A few theoretical studies of multiply 
char gecL anionic fullerenesBQ and alkali-metal (sodium) 
clusters!! have also appeared, but overall the field of mul- 
tiply anionic aggregates remains at an embryonic stage. 

In this paper, we investigate the stability and decay 
channels of Au^~ and Ag^~ at finite temperature, and 
determine their appearance sizes n 2 " (clusters with N < 
n 2 ~ are energetically unstable). Two decay channels of 
doubly anionic clusters need to be considered: (i) binary 
fission, 
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which has a welL known analog in the case of doubly 
cationic clusters JjtJ and (ii) electron autodetachment 
via emission through a Coulombic barrier Ja 
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with an analogy, to proton and alpha decay in 
atomic nuclei.E3Ej The theoretical approach we use 
is a finite-temperature semi-empirical shell-correction 
method (SCM), which incorporates triaxial shapes and 
which has been previously used successfully to describe 
the properties of neutral and cationic metal clusters .til 



Our main conclusion is-that, unlike the case of doubly 
cationic metal clusters, ElO fission of Au^ and Ag^~ is 
not a dominant process, and that the appearance sizes 
of these doubly anionic clusters are determined by elec- 
tron autodetachment. Our results for the second electron 
affinities exhibit pronounced electronic shell effects and 
are in excellent agreement with most recent experimen- 
tal datao for Au^~ with n 2 ^ — 12. For Ag^~, we predict 
r£" = 24. 

The finite-temperature multiple electron affinities of a 
cluster of N atoms of valence v (we take v = 1 for Au 
and Ag) are defined as 

A Z (N, (3) = F(/3, vN, vN + Z-1)- F(/3, vN, vN + Z) , 

(3) 

where F is the free energy, (3 = l/fegT, and Z > 1 
is the number of excess electrons in the cluster (e.g., 
the first and second affinities correspond to Z = 1 and 
Z = 2, respectively). To determine the free energy, we 
use the shell correction method. In the SCM, F is sepa- 
rated into a smooth liquid-drop-model (LDM) part -Fldm 
(varying monotonically with N), and a Strutinsky-type 
shell-correction term Af,, 



F sp - F sp , where F sp is the 
canonical (fixed TV at a given T) free energy of the valence 
electrons, treated as independent single particles mov- 
ing in an effective mean-field potential (approximated 
by a modified Nilsson hamiltonian pertaining to triax- 
ial cluster shapes), and F sp is the Strutinsky-averaged 
free energy. The smooth Fldm contains volume, sur- 
face, and curvature contributions, whose coefficients are 
determined as described in Ref. [l4|, with experimental 
values and temperature dependencies. In addition to the 
finite-temperature contribution due to the electronic en- 
tropy, the entropic contribution from thermal shaoe fluc- 
tuations is evaluated via a Boltzmann averaging.ES a ) 

We note here that the smooth contribution Az (N, /3) 
to the full muJtiple electron affinities Az(N,f3) can be 
approximated^ by the LDM expression 

~ ~ (Z-l)e 2 (Z-l + -f)e 2 

A z = At - l^ TX J - =W- y ^ /Arx V , (4) 
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where R(N) = rgN 1 / 3 is the radius of the positive back- 
ground (r s is the Wigner-Seitz radius which depends 
weakly on T due to volume dilation), 7 = 5/8, ^0 is 
an electron spillout parameter, and the work function 
W is assumed to be temperature independent [we take 
W(Au)= 5.31 eV and W(Ag)= 4.26 eV]. 
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FIG. 1. (a) Measured [see figure 3(a) in Ref. | average 
relative abundances of Au^~ clusters (i.e., the ratio of the 
number of the observed dianions over the sum of the numbers 
of corresponding singly-anionic precursors and dianions) as a 
function of cluster size. Note the logarithmic ordinate scaling, 
(b) Calculated second electron affinity (A2 in eV) for gold 
clusters at T — 300 K plotted versus N. Results from SCM 
calculations are connected by a solid line, and LDM results 
[see Eq. (^) with Z = 2] are depicted by the dashed line. 
A-2 > corresponds to stable dianionic Au^~ clusters; note 
the appearance size n a ~ = 12. Energies in units of eV. 
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FIG. 2. (a) Fission dissociation energies (Ajv,p in eV) for 
binary fission Au^~ -^Aup+Au^_ p , calculated at T = 300 
K with the SCM for parent dianionic clusters with N = 7, 14 
and 21, and plotted versus P. Note that in all cases the most 
favorable fission channel corresponds to P — 1. (b) SCM fis- 
sion dissociation energies, Ajv.i, at T = 300 K for the most 
favorable channel, plotted versus cluster size. Exothermic fis- 
sion (Ajv.i < 0) is found only for the smallest cluster. 



In a recent experiment p singly anionic gold clusters 
Au^y (N < 28) were stored in a Penning trap, size se- 
lected, and transformed into dianions, Au^T, through ir- 
radiation by an electron beam. The measured^ relative 
intensity ratios of the dianions to their monoanionic pre- 
cursors are reproduced in Fig. 1(a); they exhibit size- 
evolutionary patterns (arising from electronic shell ef- 
fects) reminiscent of those found earlier in the mass abun- 
dance spectra, ionization potentials and fipl electron 
affinities of alkali- and coinage-metal clustersJl3 Since the 
stability of the dianions relative to their monoanionic pre- 
cursors depends on the second electron affinity A%, it may 
be expected that Ai and the relative signal intensity of 
the Au^~ clusters will exhibit correlated patterns as a 
function of size. Here we note that stable dianions must 
have A2 > 0, whereas those with A2 < are unstable and 
decay via process (fU, i.e., via electron emission through 
a Coulombic barrieiu (see below). 

In He-. 1(b), we display the SCM theoretical 
resultaljO for the second electron affinity of Aujv clus- 
ters in the size range 10 < N < 30. These results cor- 
relate remarkably well with the measured relative abun- 
dance spectrum [see Fig. 1(a)]. Note in particular: (i) the 
observed and predicted appearance size n^~(Au)= 12; 
(ii) the relative instability of Au^ [portrayed by its ab- 
sence in Fig. 1(a) and the negative A 2 value in Fig. 1(b) 



associated with the closing of a spheroidal electronic sub- 
shell (containing 14 electrons) in the singly anionic Au^~ 3 
parent cluster, see Ref. |l4](b)] ; (iii) the pronounced lower 
stability of Au^g relative to its neighboring cluster sizes 
[associated with the closing of a major electronic shell 
(containing 20 electrons) in the Auj~ 9 parent cluster] ; (iv) 
the overall similarity between the trends in Fig. 1(a) and 
Fig. 1(b) (that is, even-odd alternations for N < 19 with 
a sole discrepancy at N — 15, and the monotonic behav- 
ior for N > 19). Underlying the pattern shown in Fig. 
1(b) are electronic shell effects [compare in Fig. 1(b) the 
shell-corrected results indicated by the solid line with the 
LDM curve] combined with energy-lowering shape defor- 
mations of the clusters (which are akin to Jahn- Teller 
distortions and are associated with the lifting of spectral 
degeneracies for open-shell cluster sizes). 

To explore the energetic stability of the Au^~ clus- 
ters against binary fission [see Eq. (Jjj)], we show in Fig. 
2(a) SCM results, at selected cluster sizes (N — 7, 14 
and 21), for the fission dissociation energies Ajv.p = 
F(Aup) + F(Au N _ p ) - F(Au 2 N ), with the total free en- 
ergies of the parent dianion and the singly-charged fis- 
sion products calculated at T = 300 K. For all Au^~ 
parent clusters, the energetically favorable channel (low- 
est Ajv,p) corresponds to P = 1 (i.e., one of the fission 
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FIG. 3. SCM second electron affinities [A? in (a)] and fis- 
sion dissociation energies [Ajv.i in (b)] for the most favorable 
channel (P = 1) for Ag^~ clusters at T = 300 K, plotted ver- 
sus cluster size. In (a) LDM results [see Eq. (^) with Z = 2] 
are depicted by the dashed line. Note the appearance size 
n 2_ = 24. Energies in units of eV. 

products is the closed-shell Au~ anion). The influence of 
shell effects on the fission dissociation energies is evident 
particularly in cases where the fission channel involves 
closed-shell magic products (see P = 7, and equivalently 
P = 14, for A 7 = 21, and the pronounced effect at P = 7 
for N = 14 where both fission products are magic). The 
fission results summarized in Fig. 2(b) for the most fa- 
vorable channel (P = 1) illustrate that exothermic fis- 
sion (that is Ajv,p < 0) is predicted to occur only for the 
smallest size (N — 2). This, together with the existence 
of a fission barrier, leads us to conclude that the decay 
of Au^7 clusters is dominated by the electron autode- 
tachment process (which is operative when A^< and 
involves tunneling through a Coulomb barrierO), rather 
than by fission. 

Finally, we show in Fig. 3 SCM results for the sec- 
ond electron affinity \A^ in Fig. 3(a)] and the fission dis- 
sociation energies [Ajv,p in Fig. 3(b)] corresponding to 
the most favorable channel (P = 1) for silver dianionic 
clusters Ag^~. Again binary fission is seen to be en- 
dothcrmic except for N = 2, and the appearance size for 
Ag^~ (i.e., the smallest size with A2 > 0) is predicted to 
be n^ - (Ag)= 24. The shift of the appearance size to a 
larger value than that found for gold dianionic clusters 
[that is n^(Au)= 12, see above] can be traced to the 
smaller work function of silver, as can be seen from the 
LDM curves calculated through the use of Eq. (0) with 
Z = 2.0 
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value is the same as the ratio between the experimental 
and stabilized-jellium values for the surface tension. The 
need to use experimental values for a s , a c , and W arises 
from the fact that the stabilized-jellium-LDA omits con- 
tributions from the d-electrons and thus does not provide 
accurate values for the surface tension and the work func- 
tion of gold and silver. The mass of the delocalized valence 
s-electrons was taken equal to the free-electron mass. The 
experimental temperature dependence of the surface ten- 
sion and the coefficient of linear thermal expansion were 
taken from standard tables as elaborated in Ref. |ll|(a). 
The use of a finite value [see Ref. 0(b)] for the spillout 
parameter So in the case of negatively charged clusters is 
crucial for obtaining accurate results. Taking So = (as 
was done in Ref. ^) yields a higher value for the appearance 
sizes. 



